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Abstract 

Galilei invariant equations for massless fields are obtained via con- 
tractions of relativistic wave equations. It is shown that the collection 
of non-equivalent Galilei-invariant wave equations for massless fields 
with spin equal 1 and is very broad and describes many physically 
consistent systems. In particular, there exist a huge number of such 
equations for massless fields which correspond to various contractions 
of representations of the Lorentz group to those of the Galilei one. 
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1 Introduction 



It was observed by Le Bellac and Levy-Leblond pQ already in 1973 that the 
non-relativistic limit of the Maxwell equations is not defined in a unique 
way. According to Le Bellac and Levy-Leblond ([lj p. 218), the term "non- 
relativistic" means "in agreement with the principle of Galilean relativity". 
Moreover, they claimed that there exist two Galilei invariant theories of 
electromagnetism which can be obtained by appropriate limiting procedures 
starting with the Maxwell theory. 

The combination of word "Galilean electromagnetism" itself introduced 
in pQ looked rather strange since it is pretty well known that electromagnetic 
phenomena are in perfect accordance with the Einstein relativity principle. 
However physicists are always interested when non-relativistic approxima- 
tions are adequate, which makes the results of paper [1] quite popular. The 
importance of such result is emphasized by the fact that the correct defini- 
tion of non-relativistic limit is by no means a simple problem in general and 
in the case of theories of massless fields in particular, see, for example, [2]. 

Analyzing the contents of main impact journals in theoretical and mathe- 
matical physics one finds that the interest of researches in Galilean aspects of 
electrodynamics belongs to evergreen subjects. Various approaches to Galilei 
invariant theories which were discussed briefly in [3J. An effective approach 
to Galilean electromagnetism was used in papers [I], [5] in which Maxwell 
equations were generalized in to (l+4)-dimensional Minkowski space and 
then reduced to Galilei invariant equations. Such reduction is based on the 
fact that the Galilei group is a subgroup of the generalized Poincare group 
(i.e., of the group of motions of the flat (l+4)-dimensional Minkowski space). 
For reduction of representations of the group P(l,4) to those of the Galilei 
group see refs. [B] and [7J. 

Nevertheless, the Galilean electromagnetism still contains many unsolved 
problems, for example, the question of complete description of all possible 
Galilean theories for vector and scalar massless fields. And exactly solution 
of this problem is the main issue of the present paper. 

In paper [3J indecomposable representations of homogeneous Galilei group 
HG(1,3) were derived, namely all those which when restricted to represen- 
tations of the rotation subgroup of the group HG(1,3), are decomposed to 
spin and spin 1 representations. Moreover, their connection with repre- 
sentations of the Lorentz group via the Inonii-Wigner contractions [8] were 
studied in [3J and [H] . These results open a way to complete classification of 
the wave equations describing the related fields both massive and massless 
ones. 

In the present paper we use our knowledge of indecomposable represen- 
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tations of the homogeneous Galilei algebra hg(l,3) to figure out the Galilei 
invariant equations for vector and scalar massless fields. We shall show that 
in contrast to the corresponding relativistic equations for which there are 
only two possibilities - the Maxwell equations and equations for the longi- 
tudal massless field, the number of possible Galilean equations is very huge. 
Among them there are equations with more component and less component 
fields then in the Maxwell equation. 

These results can be clearly interpreted in terms of representation and 
contraction theories. As it was proved in papers [3] and [9], there is a large 
variety of possible contractions of representations of the Lorentz group to 
those of Galilei one, and, consequently many non-equivalent Galilean mass- 
less fields. In the following sections we use these results of [3J, [S] to describe 
connections of relativistic and Galilean formulations of theories for massless 
fields. These connections appear to be rather non-trivial: in particular, com- 
pletely decoupled relativistic systems could be contracted to coupled Galilean 
ones. 

In Sections 2 and 3 we present some results of paper [3J related to the 
classification of indecomposable representations of the homogeneous Galilei 
group and contractions of the related representations of Lorentz group. These 
results are used in Sections 4 and 5 to classify all non-equivalent Galilei- 
invariant equations for massless fields. Then the classification results are 
summarized and discussed in Section 6. 

2 Indecomposable representations of the ho- 
mogeneous Galilei group 

The Galilei group G(l, 3) consists of the following transformations of tempo- 
ral and spatial variables: 

t-*t' = t + a, m 



where a, b and v are real parameters of time translation, space translations 
and pure Galilei transformations respectively, and R is a rotation matrix. 

The homogeneous Galilei group HG(1, 3) is a subgroup of G(l, 3) leaving 
invariant a point x = (0, 0, 0) at time t = 0. It is formed by space rotations 
and pure Galilei transformations, i.e., by transformations (CQ) with a = b = 0. 

The Lie algebra hg (1,3) of the homogeneous Galilei group includes six 
basis elements, i.e., three generators S a ,a = 1,2,3 of the rotation subgroup 
and three generators r\ a of Galilean boosts. These basis elements should 
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satisfy the following commutation relations 




(2) 



All indecomposable representations of HG(1, 3) which when restricted 
to the subgroup of rotations are decomposed to direct sums of vector and 
scalar representations, has been found in paper [3J. These indecomposable 
representations are labeled by triplets of numbers n, m, A and denoted as 
D(m,n, A)). The labeling numbers take the values 



In accordance with ([3D there exist ten non-equivalent indecomposable 
representations D(jn,n, A)). The spaces of representation D(m,n, A)) can 
include three types of rotational scalars A, B, C and five types of vectors 
R, U, W, K, N whose transformation laws with respect to the Galilei boost 
are: 



A -> A' = A, 

B -> B' = B + v • R, 

C^C' = C + v-XJ + \w 2 A, 

R — > R' = R, 

U -> U' = U + vA, 

W^W' = W + vxR, 

K^K' = K + vxR + vA, 

N^N' = N + vxW + vB + v(v • R) — |v 2 R 



where v is a vector whose components are parameters of the Galilei boost, 
v • R and v x R are scalar and vector products of vectors v and R. 

The carrier spaces of indecomposable representations of group hg(l, 3) 
include such sets of scalars A, B, C and vectors R, U, W, K, N which trans- 
form between themselves w.r.t. the transformations (jlj) but cannon be split 



1 < (n - m) < 2, n < 3, 
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(3) 



3 



to a direct sum of invariant subspaces. There exist exactly ten such sets: 



{^4} 
{R} 
{B,R} 

RU} 

{A,U,C} 
{W,R} 
{R,W,fl} 
{A, K, R} 
{A,B,K,R} 
{A,N,W,R} 



D(0, 
D(l, 
D(l, 
D(l, 
D(l. 



1,0 
0,0 
1,0 

1,1 
2,1 



£>(2,0,0 
■D(2,l,0 
D(2,l,l 
D(2,2,l 
D(3,l,l 



(5) 



Thus in contrary to the relativistic case where are only three Lorentz 
covariant quantities which transforms as vectors or scalars under rotations, 
i.e., the relativistic four- vector, antisymmetric tensor of second order and 
scalar, there are ten indecomposable sets of Galilean vectors and scalars 
which we enumerate in equation (jSJ). 



3 Contractions of representations of the Lorentz 
algebra 

It is well known that the Galilei algebra can be obtained from the Poincare 
algebra by a limiting procedure called "the Inonii-Wigner contraction" [8]. 
Representations of these algebras also can be connected by this kind of con- 
traction. However, this connection is more complicated in two reasons. First, 
the contraction of a non-trivial representation of the Lorentz algebra yields to 
the representation of the homogeneous Galilei algebra which the generators 
of Galilei boost are represented trivially, so to obtain a non-trivial represen- 
tation it is necessary to apply additionally a similarity transformation which 
depends on the contraction parameter in a tricky way. Second, to obtain 
an indecomposable representations of hg(l,3), it is necessary in general to 
start with completely reducible representations of the Lie algebra of Lorentz 
group. 

In paper [3] representations of the Lorentz group which can be contracted 
to representations D(m, n, A) of the Galilei group were found and the related 
contractions where specified. Here we present only a part of the results from 
[3] which will be used in what follows. 

Let us start with the representation Z?(|, |) of the Lie algebra so(l, 3) of 
the Lorentz group, whose carrier space is formed by four-vectors. Basis of 
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this representation is given by 4 x 4 matrices of the following form: 

c _ p. ( s c 3x i \ o _( 3x3 , fi x 

Oab - £abc n n , dOa ~ . n ■ W 



. o lx3 o ; ' ua v fc a 

Here s a are matrices of spin one with elements (s a )b c = ie a b c , k a are 1x3 
matrices of the form 

fci = (i,0,0), fc 2 = (0,L0), fc 3 = (0,0,i). (7) 

The Inonii-Wigner contraction consists of transformation to a new basis 

S a b > S ab} S() a >■ ES() a 

followed by a similarity transformation of all basis elements — > S' = 
yS^V" 1 with a matrix depending on contracting parameter e. Moreover, 
V should depend on e in a tricky way, such that all the transformed generators 
S' ab and eS' 0a are kept non-trivial and non-singular when e — > (BJ. 
There exist two matrices V for representation ([6]), namely 



and 



£-^3x3 O3XI 
0lx3 1 



^3X3 3 xl 



Using (jSJ) we obtain 

S' ab = VxSatVr 1 = S ab , S' 0a = e^SoaV,' 1 = ( ) . (10) 

Then, passing e to zero, we come to the following matrices 

' _ ( s a 3x i \ . „, 1 / 3x3 3x i 



S a = -e abc S bc = ^ ^ Q x j, Va = lim^J^o = Q x J • (11) 

Analogously, using matrix V2 instead of V\ we obtain 

*-(.£, v).*- )■ <-) 



Matrices (1111) and ( 1121) satisfy commutation relations p]), i.e., they realize 
representations of algebra hg(l,3). They are generators of indecopmosable 
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for £>(1,1,0) and § 



representations .D(1,1,0) and D(l, 1, 1) of the homogeneous Galilei group 
respectively. Indeed, denoting vectors from the related representation spaces 
as 

( Ri \ 

R2 
Rz 

\ B J 

and using the transformation laws (Tj0) for A, 
and U = column(L r i, U2, U3) we easily find the corresponding Galilei boost 
generators rj a in the forms (|lip and (|12p . As far as rotation generators S a are 
concerned they are direct sums of matrices of spin one (which are responsible 
for transformations of 3-vectors R and U) and zero matrices (which keep 
scalars A and B invariant). 

To obtain five-dimension representation .0(1,2,1) we are supposed to 
start with a direct sum of representations -D(|, ~) and D(0, 0) of the Lorentz 
group. The corresponding generators of algebra so(l,3) have the form 



\ A J 

B, R 



for D(l,l,l) 



column(i?i, R2, R3 



5, 







(13) 



where S^ u are matrices 



and the dots denote zero matrices of an appropri- 
ate dimension. The corresponding matrix of similarity transformation can 
be written as: 




U3 1 




(14) 



As a result we obtain the following basis elements of representation D(l, 2, 1) 
of algebra hg (1,3): 



S a 











(15) 



Matrices rj a (JTBI) generate transformations (jlj) of components of vector- 
function \1/ = column(U, A, C) so that the relation 

$ -> = exp(iT/ • v)$ 

being written componentwise, simply coincides with the corresponding trans- 
formation properties written in (|4]). 

In paper [3J was shown how all the other representations D(m,n,\) of 
homogeneous Galilei algebra can be obtained via special contractions of ap- 
propriate representations of Lorentz algebra. Here we have presented only 
such contractions which will be used later on. 
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4 Galilean massless fields 



To construct Galilean massless equations it is possible to use the same ap- 
proach which we have used in [10J to derive equations for the massive fields. 
However we prefer to apply another technique which consists in contractions 
of the relativistic wave equations. 



4.1 Galilean limits of Maxwell's equations 

According to the analysis made by Levy-Leblond and Le Bellac in 1967 [IT] . 
[JJ (see also [12], [13]) there are two Galilean limits of Maxwell's equations. 

In the so called "magnetic" Galilean limit we recieve a pre-Maxwellian 
electromagnetism. The corresponding equations for a magnetic field H and 
electric field E take the form 

VxE m -f^ = 0, V-E m = ejl 

VxH m = ej m , V-H m = 1 ; 

where j = (j m , j m ) is an electric current and e is electric charge. 

Equations ffTBT) are invariant with respect to the Galilei transformations 
(00) provided vectors H m , E m and electric current j cotransform as 

H m ► Hm; E m > E m V X I 1 ■ /-, »\ 

\ _> \ jO _> jO , • K 1 ') 

Jin jmi J m J m ' Jm- 

Introducing Galilean vector-potential A = (A , A) such as 

8A 

H m = VxA, E m = -— -VA° (18) 

at 

we conclude that the related transformation lows for A have the following 
form: 

A ^ A° + v- A, A^ A. (19) 

The other, "electric" Galilean limit of Maxwell's equations looks as 

VxH e + f = ej ej V-E e = e Je 4 , 
VxE P = 0, V • H P = 0, 



whilst the Galilean transformation lows have the following form 

H e ^H e + vxE e , E e — > E e , , - 

.4 .4 .4 l^-U 

Je -> Je + VJ e 4 , J 4 J 4 . 
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Vectors H e and E e can be expressed via vector-potentials as follows 

H e = V x A, E e = -VA 4 (22) 



and so the corresponding Galilei transformations for the vector-potential are: 
A 4 -> A\ A -> A + vA 4 ; (23) 

Till this point we just presented the results of the Levy-Leblond anal- 
ysis of two possible Galilean limits of Maxwell's equations. Let us note 
that these results admit clear interpretation using the representation theory. 
Indeed, in accordance with the results presented in Section 3 there exist ex- 
actly two non-equivalent representations of the homogeneous Galilei group 
the carrier spaces of which are four- vectors - the representations Di(l, 1, 0) 
and Di(l, 1, 1). In other words there are exactly two non-equivalent Galilei 
transformations for four-vector-potentials and currents, which are given ex- 
plicitly by equations ( 1T91) . (TT71) and ( 1231) . ( 12TT) . Equations for massless fields 
invariant with respect to these transformations are given by relations ( 1T61) 
and f[2"U]) respectively. 

Both representations, i.e., Di(l, 1, 0) and Di(l, 1, 1), can be obtained via 
contractions of the representation D(l/2, 1/2) of the Lorentz group whose 
carrier space is formed by relativistic four-vectors. The related contraction 
matrices are given explicitly by relations (E]) and ([9]). Each of these contrac- 
tions generates the Galilean limit of Maxwell's equations, and in this way 
we obtain easily the systems (I16p and (1201) . In the following section we will 
obtain equations (TTtjj) and (I2T)]) via contraction of a more general system of 
relativistic equations for massless fields. 

4.2 Extended Galilean electromagnetism 

In accordance with our analysis of vector field representations of the Galilei 
group there exists the only representation, namely, _Di(l,2, 1) whose carrier 
space is formed by five-vectors. Such five-vectors appears naturally in many 
Galilean models, especially in those ones which are constructed via reduction 
approach |3j, i.e., starting with models invariant with respect to the extended 
Poincare group P(l,4) and making reduction to its Galilean subgroup. 

As mentioned in [3], there is a formal possibility to introduce such five 
component gauge fields which join and extend the magnetic and electric 
Galilean limits of relativistic four-vector potentials. However the physical 
meanings of the related theories was not clear. Moreover, Maxwell's elec- 
trodynamics can be contracted either to the magnetic (fTol) or to the electric 
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limit fl2"U]) . and it is generally accepted to think that it is impossible to for- 
mulate a consistent theory including both 'electric' and 'magnetic' types of 
Galilean gauge fields, see, e.g., [5]. 

In contrary to [5], we shall show that it is possible to join the 'electric' 
and 'magnetic' Galilean gauge fields since the Galilean five-vector potential 
appears naturally via contraction of a relativistic theory. Rather surprisingly, 
the related relativistic equations can be decoupled to two non-interacting 
subsystems whereas its contracted counterpart appears to be coupled. This is 
in accordance with the observation presented in [3] that some indecomposable 
representations of the homogeneous Galilei group appears via contractions 
of the completely reducible representations of the Lorentz group. 

Let us start with relativistic equations for the vector-potential 

p^A" = -ef (24) 

in the Lorentz gauge 

p^A* = or p A° = p A. (25) 

Consider also the inhomogeneous d'Alembert equation for a relativistic 
scalar field which we denote as A 4 : 

p^p.A* = ef. (26) 

Introducing the related vectors of the field strengthes in the standard 
form 

H = V x A, E = -ip - VA°, F = VA\ F° = ^- (27) 

ox ox 

we come to Maxwell's equations for E and H: 

OH 

V x E - — = 0, V - H = 0, 

dx 



:28) 



d~E 

VxH+ — = ej, V-E = ej° 

dx 

and the following equations for F and F° 
OF 



OX: 



+ V • F = ef, 



o 



OF 

V x F = 0, — = VF°. 

OXn 



(29) 
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Surely the system of equations (jzoj) and (|29|) is completely decoupled. 
Rather surprisingly its Galilean counterpart which will be obtained using 
the Inonii-Wigner contraction appears to be coupled. This contraction can 
be made directly for equations ( 1281) . ( 1291) but we prefer a more simple way 
with using the potential equations ( |24|) . 

The system of equations (!!M]) - (p3]) is a decoupled system of relativistic 
equations for the five component function 

A = column^ 1 , A 2 , A 3 , A , A 4 ) = column(A, A , A 4 ). (30) 

Moreover, the components (A 1 , A 2 , A 3 , A ) transform as four- vector and A 4 
transforms as a scalar, so the generators of the Lorentz group defined 
on A are direct sums of the four-vector generators S^ v and zero matrices 
given by equation ( fl3|) . In other words, such generators realize the direct 
sum -D(|, |) © D(0, 0) of representations of algebra so(l, 3). 

In accordance with results of paper [3] the completely reduced represen- 
tation of the Lie algebra of Lorentz group whose basis elements are given 
by equation (fl3l) can be contracted either to direct sum of indecomposable 
representations of the the Galilei algebra hg(l, 3) or to indecomposable rep- 
resentation Di(l,2, 1) this algebra. 

Let us consider the first possibility, i.e., the contraction to the indecom- 
posable representation. Such contraction is presented in Section 3, see equa- 
tions (fl3l-(TT5!) here. 

Let us demonstrate that this contraction reduces the decoupled relativis- 
tic system ( |24|) . (|26|) to a system of coupled equations invariant with re- 
spect to the Galilei group. Indeed denoting A 1 = U 3 A = (A', A' , A' A ) and 
j' = U^j = j' 4 ) and taking into account that the relativistic variable Xq 
is related to non-relativistic variable t as xq = ct, we come to the following 
system of equations for the transformed quantities: 

P 2 A' k = -ef, i^ = p-A'. (31) 

Generators of Galilei group for vectors A' and j' are given by equation 
(I15p . so under the Galilei transformation ^TJ) they cotransform in accordance 
with the representation Dx(l,2, 1), i.e., 

A -> A + v • A + —A 4 , A -> A + vA 4 , A 4 -> A 4 , (32) 



and 



2 



J 4 - /, J - j + vj 4 , f - j° + v • j + ^ 2 j 4 . (33) 
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Of course transformations (pQ), (jH2D and f[3"3l keep the system f[3~Tl) invari- 
ant. In accordance with f[3"Tj) the related field strengthes (compare with (12"T|) ) 



W = V x A', N 



<9A' 



/0 



R = VA 



n 



B 



dA 



(34) 



dt 



dt 



satisfy the following equations 



J° = V • N - gB - ej° = 0, 
J=VxW + VB-ej = 0, 
J 4 = V ■ R - ej 4 = 0, 
Af = |W + V x N = 0, 
W = |R-V5 = 0, 
UeVxR=0, 
£ = V- W = 0. 



(35) 



Like (EH), equations (1351) are covariant with respect to the Galilei group. 
Moreover, the Galilei transformations for fields R, W, N , B and current 
j are given by equations (T4|) and fl33l respectively. In other words, these 
fields and the current j form carrier spaces of representation -Di(3, 1, 1) and 
Di (1,2,1) of algebra hg (1,3) correspondingly. 

In contrast with the decoupled relativistic system of equations ( 1281) and 
( |29l) its Galilean counterpart ( 1351) appears to be a coupled system of equations 
for vectors R, W, N and scalar B. 

The system of equations (!33|) was obtained in paper [1] by reduction of 
generalized Maxwell equations invariant with respect to extended Poincare 
group P(l, 4) including one time variable and four spatial variables. We prove 
that this system is nothing but the contracted version of system ( |28l) . ( 1291) 
including the ordinary Maxwell equations and equations for four-gradient of 
scalar potential. It other words, the system of Galilei invariant equations 
(133|) admits a clear physical interpretation as a non-relativistic limit of the 
system of familiar equations ( |28l) and ( !29l) . 

4.3 Reduced Galilean electromagnetism 

In contrast with the relativistic case the Galilei invariant approach makes it 
possible to reduce the number of field variables. For example, considering the 
magnetic limit (fl6l) of the Maxwell equations it is possible to restrict ourselves 
to the case H m = in as much as this condition is invariant with respect to 
the Galilei transformations ffTT|) . Notice that in the relativistic theory such 
condition can be imposed only in a particular frame of references and will be 
affected by the Lorentz transformation. 
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In the mentioned sense equations flHoj) are reducible too. They are defined 
on the most extended multiplet of vector and spinor fields which is a carrier 
space of indecomposable representation of the homogeneous Galilei group. 
The related representation -Di(3, 1, 1) is indecomposable but reducible, i.e., 
includes subspaces invariant with respect to the Galiei group. This makes 
it possible to reduce the number of dependent components of equations (J35J) 
without violating its Galilei invariance. 

In this section we consider systematically all possible Galilei invariant 
constrains which can be imposed on solutions of equations ( 1351) and present 
the related reduced versions of Galilean electromagnetism. 

In accordance with (TjJ and fl33|) the vector R and the fourth component j 4 
of current form invariant subspaces with respect to Galilei transformations. 
Thus we can impose the Galilei-invariant conditions 

R = or VA 4 = 0, f = (36) 

and reduce the system (135]) to the following one 



P + VxE = 0, 

V x H = ej, V • H = 0, , , 

V-E = |5 + ej°, {67} 

vs = o 

where we denote H = W|r= , E = N|r= , 5* = -£>|r=o- 

The vectors E, H and scalar S belong to a carrier space of the represen- 
tation Di(2, 1, 1). Their Galilei transformation laws look as 

E^E + vxH + vS, 

H-H, S^S. (38) 

In accordance with ( |38i) S belongs to invariant subspace of Galilean trans- 
formations, so we can impose one more Galilei invariant condition 

S = or -— = 0. (39) 
dt v ; 

As a result we come to equations (fT6l) . i.e., to the magnetic limit of Maxwell's 
equations. Thus equations (ITB"]) are nothing but the system of equations (l3"5"j) 
with the additional Galilei invariant constrains (1361) and (1391) . 

Galilei transformations for solutions of equations (Tl6l) are given by equa- 
tions (fl71) . Again we recognize the invariant subspace spanned on vectors 
H m , and so it is possible to impose the invariant condition 

H m = or A = Vip, j m = (40) 
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where (p is a solution of the Laplace equation. As a result we come to the 
following system 

V x E = 0, V • E = ep (41) 

where we denote E = E m | Hm = and p = j° m \u m =o- 

Equation (j4~Ii) is still Galilei invariant, moreover, both E and p are not 
changing under the Galilei transformations. The related potential A is con- 
strained by conditions (|36|) . (|39|) and (|40|) . Moreover, up to gauge transfor- 
mations it is possible to set in (1401) A = 0. 

4.4 Other reductions 

Equations (1371) . (fl6l) and (f4TT) exhaust all Galilei invariant systems which can 
be obtained starting with (IHoT) and imposing additional constraints which 
reduce the number of dependent variables. To find other Galilei invariant 
equations for massless vector fields we use the observation that the Galilei 
transformations for equations ( 135]) have the form given by relations (T41) and 
( |33l) if we change N — > A/", W — > W, • • • here. Thus the subsystem of 
equations (1331) obtained by excluding equations Af = and J7" = 0: 

J = V x W + f R - ej = 0, 
J 4 = V ■ R - ej 4 = 0, 

W = |R-Vfi = 0, (42) 
^ = V x R = 0, 
£ = V-W = 

is Galilei covariant too and does not include dependent variables N and 
j°. The Galilean transformations for W, R, B and j,j 4 are still given by 
equations (j4j) and (1331) . 

Following the analogous reasonings it is possible to exclude from fj4~2"]) 
equations J~ = and £> = and obtain the system 

V ■ R - ej 4 = 0, 

fR-V£? = 0, (43) 

V x R = 

which includes only two vector and two scalar variables. The related potential 
without loss of generality reduces to the only variable A 4 . 

The other way to reduce system (T42l) is to exclude the equation W = 0. 
As a result we come to the electric limit for the Maxwell equation (1201) for 
W = H e and R = E e . 
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Thus in addition to flHol) . flH7|) . ffTol) and (jHJ) we have three more Galilei 
invariant systems given by equations fp£Zl) . f|4*H|) and (12U1) . These equations 
admit additional reductions by imposing Galilei invariant constrains to their 
solutions. 

Considering ( |20l) we easily find a possible invariant conditions E e = 
0, j 4 = which reduce it to the following equations 

V x H = ej, V • H = (44) 

where we denote H = H e | Ei .= . 

Let us return to system (14"2I) . This system can be reduced by imposing 
the Galilei invariant conditions R = 0, j 4 = to the following form: 

Vx *- ej = °' (45) 
V-H = 0, V5 = K ' 

where we change the notations W — > H and B — ^ S. The related potential 
have the form A = (A 4 ,0, A), moreover, A 4 should satisfy the condition 
VA 4 = 0. 

We see that in contrast with the relativistic theory there exist a big variety 
of equations for massless vector fields invariant with respect to the Galilei 
group. The list of such equations is given by formulae ffTol) . fEUl) . flH5"j) . flHTj) . 



5 Discussion 

The revision of classical results [i] related to Galilean electromagnetism made 
in the present paper appears to be successful. Our knowledge of indecom- 
posable representations of homogeneous Galilei group defined in vector and 
scalar fields [3j made it possible to complete the results of Le Bellac and 
Levy-Leblond [1] and present an extended class of Galilei-invariant equations 
for massless vector fields. Among them are decoupled systems of first order 
equations which include the same number of components as the Maxwell 
equations and equations with other numbers of components as well. The 
most extended system includes ten components while the most reduced one 
has three components only. 

It is necessary to stress that the majority of obtained equations admit 
clear physical interpretations. Thus equation (T4T]) and (jHJ) are applied in 
electro- and magnetostatics respectively. The very procedure of deducing of 
Galilei invariant equations for vector fields used in the present paper makes 
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their interpretation to be rather straitforward, since any obtained equation 
has its relativistic counterpart. 

We see that the number of Galilean wave equations for massless vector 
fields is rather extended, and so there are many possibilities to describe 
interaction of non-relativistic charged particles with external gauge fields. 
Some of these possibilities have been discussed in [3] and [TU], se also [7j, [H] 
and [I], [SJ. Starting with the found equations and using the list of functional 
invariants for Galilean vector fields presented in [9] it is easy to construct 
nonlinear models invariant with respect to the Galilei group, including its 
supersymmetric extensions. 

Invariance of any particular found equation with respect to the Galilei 
transformations can be verified by direct calculation. The main result pre- 
sented in this paper is the completed description of all such equations. 
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